Phase diagram of a strongly interacting polarized Fermi gas in one dimension 
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Based on the integrable Gaudin model and local density approximation, we discuss the ground 
state of a one-dimensional trapped Fermi gas with imbalanced spin population, for an arbitrary 
attractive interaction. A phase separation state, with a polarized superfluid core immersed in an 
unpolarized superfluid shell, emerges below a critical spin polarization. Above it, coexistence of 
polarized superfluid matter and a fully polarized normal gas is favored. These two exotic states 
could be realized experimentally in highly elongated atomic traps, and diagnosed by measuring 
the lowest density compressional mode. We identify the polarized superfluid as having an FFLO 
structure, and predict the resulting mode frequency as a function of the spin polarization. 

PACS numbers: 03.75.Ss, 05.30.Fk, 71.10.Pm, 74.20.Fg 



Strongly attractive Fermi gases with imbalanced 
spin components are common in different branches of 
physics [ll, so spin-polarized ultra- 

cold atomic gases are an atomic analog of many other 
exotic forms of matter. The true ground state of an 
attractive polarized Fermi gas remains elusive, because 
the standard Bardeen-Cooper-Schrieffer (BCS) mecha- 
nism requires the pairing of two fermions with opposite 
spin. Polarized Fermi gases cannot be explained by the 
BCS theory, due to mismatched Fermi surfaces. Various 
exotic forms of pairing have been proposed, including 
the deformed Fermi surface [3] , interior gap pairing [1] or 
Sarma superfluidity [H] , phase separation [(| , and the in- 
homogeneous Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state [3| . 

Recent measurements on polarized 6 Li gases [HQ near 
a Feshbach resonance provide a route towards testing 
these theories in experiment. However, the presence of 
a harmonic trap in these experiments makes it difficult 
to identify which if any of the different pairing schemes 
occurs. Theoretically, it is desirable to have an exactly 
soluble mode of polarized uniform Fermi gases to identify 
various pairing schemes, and clarify the issue of the trap. 

In the present Letter we report on the exact ground 
state of a homogeneous ID polarized Fermi gas 
with attractive inter-component interactions at zero tem- 
perature. We then study the phase diagram of an inho- 
mogeneous Fermi gas under harmonic confinement, using 
the local density approximation (LDA). We complement 
this with a mean-field Bogoliubov-de Gennes (BdG) the- 
ory in the weak-coupling limit, where the phase fluctua- 
tions are small, in order to clarify the physical meaning 
of the solutions. Collective mode frequencies are also 
calculated as an experimental diagnostic, thus extending 
previous work on the unpolarized case 

As well as being a theoretical test bed for the ground- 
state issue, ID Fermi gases in traps can be realized using 
two-dimensional optical lattices, where the radial motion 



of atoms is frozen out. Thus, one can experimentally 
check these quantum many-body predictions, which has 
also been recently carried out for the ID Bose gas flj| . 
The following remarkable features are found: 

(A) In the ground state of a uniform system, we find 
three distinct phases with increasing chemical potential 
difference between species: an unpolarized BCS super- 
fluid, a polarized superfluid, and a fully polarized normal 
state. The polarized superfluid is most widespread, and 
reduces to the FFLO-type for weak interactions. There- 
fore, it is relatively easy to observe the FFLO physics in 
ID, as anticipated in previous approximate studies fl(i |. 
In earlier work the phase diagram was not conclusive, as 
the nature of the transition from BCS to FFLO states 
was under debate 

(B) Within the local density approximation (LDA), 
we consider the inhomogeneous phase diagram of the 
trapped gas. This leads to a phase separation, as the 
inhomogeneous cloud separates into either a mixture of 
a polarized superfluid and an unpolarized superfluid, or 
a polarized superfluid and a fully polarized normal gas. 

(C) We calculate the longitudinal size of the two spin 
components and the frequency of the lowest density com- 
pressional mode. These quantities give a measurable fin- 
gerprint of the whole phase diagram. 

We describe a ID polarized Fermi gas with N = iVf + 
N± fermions each of mass m and spin polarization P = 
(iVf — Ni)/N > in a harmonic trap, by 
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is the Hamiltonian of Gaudin model of a spin 1/2 
Fermi gas attracting via a short range potential gioS(x) 
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Figure 1: (Color online) The energy per particle, chemical po- 
tential for spin up atoms (a), and chemical potential difference 
(b) in units of h 2 n 2 /2m as a function of the spin polarization 
at two coupling parameters 7 = 0.5 and 7 = 2.0. 

[j| . The coupling constant gm (< 0) can be ex- 
pressed through the ID scattering length am, gm = 
— 2h 2 / (maiu). A two-fermion bound state arises once 
N[ > 0, with binding energy £b = h 2 /(ma 2 D ). 

In the absence of the harmonic trap, the integrable 
Gaudin model, Eq. {2]), can be solved exactly using 
Bethe's ansatz 0, [lfj. Introducing linear number den- 
sities, n = N/L and n a = N a /L (a =1,4), where L is 
the size of the system, the uniform gas is characterized 
by the polarization p = (rij — ri\)/n > and a dimen- 
sionless parameter 7 = -mjic/(?i 2 n) = 2/ (nam) > 0. 
The ground state is obtained by numerically solving the 
Gaudin integral equations @, We have clarified the 
physical nature of the resulting solutions by also solving 
the weak-coupling mean-field BdG equations. 

Fig. 1 shows the energy per particle, Eh om /N, 
the chemical potential of spin up fermions, /J.hom,T = 
dEbom/dNj, and the chemical potential difference, 
5/J-hom — dEi lonl /d(N^ — Ni), as a function of polariza- 
tion p at two intermediate interaction strengths 7 = 0.5 
and 7 = 2.0. The mean-field calculations lead to the 
same general behavior. Their asymptotic behavior in the 
weak and strong coupling limits may be understood ana- 
lytically. Of particular interest is the chemical potential 
difference. In the weakly interacting limit of 7 <C 1, 
Sfihom is given by (7 < max{p, 1 - p}) 



tin 7r ti n 
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(3) 



where the first term on the right hand side is the Fermi 
energy difference of an ideal polarized gas, and the sec- 
ond term arises from the mean-field Hartree-Fock inter- 
actions. The chemical potential difference increases with 
increasing 7, and reaches a half of the binding energy of 
bound states in the strongly attracting limit of 7 — * 00, 
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Figure 2: (Color online) Uniform phase diagram, displaying 
"normal" (N), unpolarized superfluid (SF) and polarized su- 
perfluid states (SFp). At weak coupling or high density, the 
predominant SFp phase corresponds to the mean-field FFLO 
solution . The dashed and dotted lines are asymptotic behav- 
ior as described in the text. 



The first two terms in magnitude coincide with the chem- 
ical potential of a Tonks-Girardeau bosonic gas of paired 
Ni dimers flil . QjJ], which is fermionized due to strong 
attractions, while the third term is equal to the chemical 
potential of residual unpaired iVj — iVj. non-interacting 
fermions. Therefore, in the strong coupling regime the 
system behaves like a coherent mixture of a molecular 
Bose gas and fully polarized single-species Fermi sea. 

We analyze the phase structure of the ground state of 
uniform polarized Fermi gases (Fig. 2). Given an inter- 
action strength the chemical potential difference at finite 
polarization p is bounded by two critical values, S^ c ,p=o 
and Sfx c , P =i, as indicated by arrows in Fig. lb for 7 = 2.0. 
Below <fyt CiP= o, the gas remains in the BCS-like superfluid 
state with zero polarization (SF), while above 5(i c . p= i, a 
fully polarized normal state is favored (N). In between, 
a superfluid state with finite polarization (SFp) domi- 
nates. The mean-field calculation of a uniform gas shows 
that the SFp is of FFLO character, as the exactly-soluble 
ground state energy corresponds precisely with the FFLO 
solution in the weak-coupling limit. Physically S^ c ,p=o is 
the energy cost required to break pairs in unpolarized 
superfluid, i.e., the spin gap [HI, while £// c ,p=i is also 
associated with the pair-breaking (for the last pair), but 
is promoted upwards due to the Pauli repulsion from ex- 
isting fermions. The dependence of S^ c ,p=o and <5/^ CjP= i 
on the parameter 7 is reported in Fig. 2, resulting in a 
homogeneous phase diagram. 

In units of H 2 n 2 /2m, these have the following limiting 
behavior: <5/x C:P= o ~ 2y/iFf exp[— ir 2 /2 r y] and 5^ c ,p=i — 
(ir 2 1 2) as 7 — * 0, while £/i c ,p=o — e &/2 + 7r 2 /16 and 
<fy*c,f)=i — e&/2 + 7r 2 as 7 — * 00. Both critical chemical 
potential differences saturate to the half of the binding 
energy in the strong coupling limit. Converting 7 into the 
chemical potential /ihom = dEi lom /dN , we can obtain the 
phase diagram in the /Xhom-^Mhom plane. 
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Figure 3: (Color online) Density profiles of each spin com- 
ponent and their difference at a typical interaction coupling 
of Na\ D /a 2 ho = 1 for the SFp-SF phase (a) and the SFp-N 
phase (b). 

We now turn to describe the ID polarized gas in a har- 
monic trap. We partition the system into cells that can 
be treated locally as being uniform, with a local chem- 
ical potential. This local density approximation (LDA) 
is applicable provided that the number of fermions in a 
cell is much greater than unitary, and the variation of 
the trap potential across the cell is small compared with 
the local Fermi energy and hence the interface effects are 
negligible [l3|. Overall it requires a sufficient large local 
density, which implying N >■ 1, a condition readily sat- 
isfied in the ID experiment. Note that the breakdown 
of LDA has been observed in the elongated 3D trap [l7| . 
when the linear density in the transverse axis becomes 
small. The LDA amounts to determining the chemical 
potential fx = (/if + /ij.)/2 and the chemical potential dif- 
ference Sfi = (/ij — /ix)/2 of the inhomogeneous gas from 
the local equilibrium condition, 

Mhom,cr [n(x),p(x)} + \^nuj\ Q x 2 = [i a . (5) 

The normalization conditions are N ~ J n(x)dx and 
NP = J n(x)p(x)dx, where n(x) is the total linear num- 
ber density and p(x) the local spin polarization. By 
rescaling the chemical potentials, coordinate and lin- 
ear density into dimensionless form [jl a — [i a jtb, x = 
x /( a ho/ a iD) and n = nauy], the normalization equa- 
tions can be rewritten as Na^/a^ — j h(x)dx and 
i.Na\ D /a 2 ho )P = J n(x)p(x)dx. These expressions em- 
phasize that the dimensionless coupling constant in a 
trap is controlled by Na\ D /a\ , where Na 2 D /a 2 a » 1 
corresponds to weak coupling while Naf D /af lo <C 1 cor- 
responds to the strongly interacting regime. 

The qualitative feature of density profiles n a (x) is sim- 
ple to understand. Within the LDA, the local chemical 
potential n(x) decreases parabolically away from the cen- 
ter of the trap while the local chemical potential differ- 
ence 5[i(x) stays constant. It is then clear from Fig. 2 
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Figure 4: (Color online) (a) Phase diagram of a inhomoge- 
neous ID polarized Fermi gas. The dashed and dotted lines 
are asymptotic limits described in the text, (b) The Thomas- 
Fermi radius of each spin component as a function of the 
polarization at Na\ D /a\ a = 0.01, 1, and 100. 



that with a nonzero P we always have a polarized su- 
perfluid in the center where the local chemical potential 
(interaction parameter) is large (small). Away from the 
center with decreasing local chemical potential, the Fermi 
gas goes into either an unpolarized superfluid (<5/i < £b/2) 
or a fully polarized normal cloud (6[i > £b/2). Thus, 
there is a critical chemical potential difference 8^ c = et/2 
that separates the inhomogeneous system into two phase 
separation states: a mixture of a polarized superfluid core 
and an unpolarized superfluid shell (SFp-SF), or a coex- 
istence of a polarized superfluid at the center and a fully 
polarized normal gas outside (SFp-N). 

The former phase is exotic, as the BCS-like superfluid 
state occurs at the edge of the trap, in marked contrast 
to the 3D case [H]. This is caused by the peculiar effect 
of low dimensionality, for which the gas becomes more 
nonideal with decreasing ID density towards the edge of 
the trap, and hence the energy required to break the pairs 
approaches efc/2 from below. As <5/i < eh/2, there should 
be a fully paired region once the local critical chemical 
potential 5fi c , P =o > 6fi, i.e., the BCS-like superfluid. 

We show the density profile of each component in Fig. 
3 with a typical coupling parameter Nof-p/a^ — 1. In 
addition, we have performed a BdG calculation for the 
trapped gas and observe a FFLO-N phase. The resulting 
density profiles on the weak coupling side are in perfect 
agreement with the LDA calculation, indicating unam- 
biguously that the SFp phase is a FFLO state. 

We determine the phase diagram of the inhomogeneous 
polarized ID Fermi g clS clS el function of the interaction 
strength and spin polarization by calculating the criti- 
cal polarization P c that corresponds to <5/z c as a func- 
tion of the coupling constant, and plot it in Fig. 4a. 
The asymptotic behavior of P c can be computed analyt- 
ically in the weak and strong coupling limits. We find 
that P c ~ ^/{Na\ D /a 2 ho ) as Na\ D /a 2 ho -> oo, and P c ~ 
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Figure 5: (Color online) Square of the lowest density breath- 
ing mode frequency, uj 2 , as a function of the spin polarization 
for interaction parameters as indicated. The inset shows the 
mixing angle that minimizes the sum-rule mode frequency at 
Na 2 1D /a 2 ho = 1.0. 



1/5 - (256/225^X0.4^/^)1/2 as Na 2 D /a 2 ho - 0. 

We consider the experimental relevance of the two 
phase separation states, by calculating the size of the 
cloud and the lowest density compressional mode. These 
are readily detectable via absorption imaging. Fig. 4b 
reports the evolution of the Thomas-Fermi radius of two 
spin components as a function of polarization at three 
different interaction couplings. The radius for spin up 
and down fermions is the same in the SFp-SF phase, but 
diverges in opposite directions in the SFp-N phase. 

Using a sum rule approach, the frequency ui of the low- 
est density and spin density compressional (breathing) 
modes of ID trapped gases can be calculated from the 
identity h 2 uj 2 = mi/m_i, where mi = {[F+, [H,F]]) /2 
and m_i = — x(F)/2 are two energy-weighted moments 
of the breathing operator F — cos 9 ^ x 2 ^ +sin 9 J^j x ji > 
with a mixing angle — ir/2 < 9 < ir/2. The linear static 
response of the system \ can be calculated in terms 
of the susceptibility matrix ! £hom/ 'dn a dn a ' . Basically 
these two different breathing modes correspond to the 
in-phase and out-of-phase oscillations of the two imbal- 
anced spin populations. They are decoupled for a mixing 
angle 6i„ > or 9 out < that minimizes the mode fre- 
quency, analogous to the spin-charge separation in ID. 
Then, the operators F with these angles are anticipated 
to exhaust all the weights in the density and spin density 
channels, respectively. Thus, the sum rule approach is 
well applicable, providing only an upper bound on the 
mode frequency. A stringent test of this single mode 
approximation merits further study, e.g., by using the 
random-phase approximation theory. Fig. 5 shows the 
density breathing mode frequency as a function of polar- 
ization at several interaction strengths. With increasing 
polarization the frequency initially rises in the SFp-SF 
state and then gradually decreases to the ideal gas re- 
sult 2u>ho in the SFp-N phase. A peak structure is found 
that gives an independent means of identifying the FFLO 



phase, which dominates the phase diagram in the vicinity 
of this peak. 

We emphasize that the FFLO superfluid reported here 
can be detected via a Josephson junction that is formed 
by confining two ID polarized gases in a double- well trap. 
There is also a signature present in the density correla- 
tion function, which has an oscillatory behaviour in this 
phase. Further details will be provided elsewhere. 

In conclusion, we have investigated a ID polarized 
Fermi gas in a harmonic trap, and have shown that the 
trap generally gives rise to phase separation: with at 
least one FFLO-type phase present. Two distinct exotic 
phase-separated structures can occur, and are detectable 
via absorption imaging and collective mode experiments. 
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search Council Center of Excellence and by the National 
Science Foundation of China under Grant No. NSFC- 
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Note added. — After this work was completed, we be- 
came aware of a related work of G. Orso |18|. 
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